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Here, we address the question of the validity of an effective description for hyperbolic metamate- 
rials in the near- field region. We show that the presence of localized modes such as surface waves 
drastically limits the validity of the effective description and requires revisiting the concept of ho- 
mogenization in near- field. We demonstrate from exact calculations that one can find surface modes 
in spectral regions where the effective approach predicts hyperbolic modes only. Hence, the pres- 
ence of surface modes which are not accounted for in the effective description can lead to physical 
misinterpretations in the description of hyperbolic materials and their related properties. 
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With today's nanotechnology it is possible to manu- 
facture artificial composite materials at tiny scales which 
manifest optical properties that we do not encounter in 
nature. These media also called metamaterials are usu- 
ally structured at the length scale or below the wave 
length of photons. These media have been largely ex- 
ploited during the last decades to overthrow long stand- 
ing paradigms in physics. For instance, new phenomena 
such as the negative refraction [l[ , the super-resolution Q 
or reversed Doppler effects ^ were predicted for meta- 
materials. 

One class of metamaterials, the so called indefinite or 
hyperbolic materials (KM) has recently attracted much 
attention. In fact, the dispersion relation of electromag- 
netic waves in these media can be represented by hy- 
perbolic isofrequency curves meaning that they support 
propagating modes having extremely large wavevectors 
far beyond the light line [4, 5]. Because of that property 
HMs are viewed as promising candidates to develop or 
improve breakthrough technologies among which are the 
sub-diffractive imagery (hyperlensing) [6|-l8|, the thermal 
management at nanoscale (9l-[ll|, the near-field energy 
12 1, enhanced light emission [l3|, and quan- 



neous emission it was shown that the EMT can greatly 
underestimate or overestimate 16-181 the Pur cell factor 



conversion 



tum information systems [l4J. 

The exotic properties of HMs in the near-field regime 
are conventionally investigated within the framework of 
the effective medium theory (EMT). Naturally the ques- 
tion arises under which conditions the EMT remains valid 
in the close vicinity of a composite material. In general it 
is believed that for propagating modes the EMT is valid 
as long as the wavelength is much larger than the size of 
the unit-cell of the composite material. In the evanescent 
regime this condition imposes that the wavevector of the 
evanescent field has to be smaller than the inverse of the 
unit-cell size as well [15]. In the context of the sponta- 



when this condition is not met. In addition, when the 
composite materials support surface waves one encoun- 
ters a nonlocality which cannot be adequately described 
by the EMT [11, [TH. Another nonlocal effect can arise 
from the nonlocal response of the composite materials 
for very large wave vectors [2o|. Therefore, a nonlocal 
effective description is needed to take such effects into 
account 



2ii. 



In this letter, we show that the presence of sur- 
face waves at the interface of the basic components of 
HMs also requires revisiting the concept of homoge- 
nization which otherwise, in the EMT framework, can 
lead to quantitatively wrong and physically misinter- 
preted results if one estimates the hyperbolic bands from 
EMT [22, J3j . By comparing the exact calculation of the 
local density of states of electromagnetic field at arbi- 
trary distances above an HM with the predictions given 
by the EMT we highlight the regions where the effec- 
tive theory fails to describe properly the field outside the 
HM. In particular, we demonstrate from exact calcula- 
tions the existence of surface modes inside the hyperbolic 
bands which are not predicted at all by the EMT. Hence, 
some properties attributed to the hyperbolicity of mate- 
rial could be misinterpreted or even worse could not exist 
anymore. For example large changes in the Purcell fac- 
tors within the hyperbolic frequency bands could be due 
to such surface modes, which are not taken into account 
in the EMT. Hence, the EMT would in this case give 
results which are quantitatively and qualitatively wrong. 

To estimate the range of validity for the EMT in the 
presence of surface modes we consider a ID periodic 
structure composed by two materials with one which nat- 
urally supports surface waves in the spectral range of in- 
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Figure 1: (Color online) (a) Sketch of the ID SiC/Ge periodic 
structure with their layer thicknesses h and h. The period 
of the structure is a — h -\- h and the number of layers is 
equal to = 100. The surrounding medium is vacuum, (b) 
Real part of the effective permittivities e|| and e± for / = 0.4. 
The grey areas Ai and A2 mark the hyperbolic regions where 
Re(e||)Re(e±) < 0. 



terest. This structure is depicted in Fig. [TJa). One basis 
material is Silicon carbide (SiC) and the other one is ger- 
manium (Ge). The polar medium SiC supports surface 
phonon polariton (SPP) modes in the mid-infrared at 
Aspp = 10.3 fim and its permittivity is given by 



(1) 



following the Drude Lorenz model [2^ with Coo = 6.7, 
CJL = 182.7 • lO^^rad/s, ujt = 149.5 • lO^^^ad/s and 
7 = 0.9 • lO^^rad/s. The dielectric permittivity of Ge 
in the same frequency range is constant and equal to 
ece = 16. According to the EMT, for a < Aspp the 
equivalent homogenized medium to this structure is an 
uniaxial anisotropic medium with an effective permittiv- 
ity parallel and per2endicular to the optical axis given 
by the expressions [25 1 



^11 



/esiC + (l-/)eGe and e± 



esicece 



fcGe + (1 - /)eSiC 

(2) 

with the volume filling fraction / = li/a. The reflection 
coefficients for s- and p-polarized waves in the effective 
description are thus 



kzo 
kzo + 



and 



kn 



e\\kzo + kp 



(3) 



with 

kr, = 



HI 



.2^ 



K?^ ks = Y^cj^/c^eii — tz'^ and 
the normal components of the 



wavevector in vacuum and in the medium for s and p 
waves. Here the wave vector = {kx^kyY is perpendic- 
ular to the optical axis. Therefore p-polarized electro- 
magnetic waves (the so called extra-ordinary waves) in 
such an effective uniaxial medium fulfill the dispersion 
relation [26| 



to 



r.2 ■ 



(4) 



If Re(e||) and Re(e^) are both positive Eq. (|4]) describes 
an elliptic dispersion curve in the {n^kz) plane while on 
the contrary if the effective permittivities have opposite 
signs, i.e., if Re(e||)Re(e^) < 0, the dispersion curve de- 
scribed by Eq. (H]) becomes hyperbolic. In the spectral 
range where Re(e||)Re(e^) < waves are propagating 
and are called hyperbolic modes. In Fig.[TJb) Re(e||) and 
Re(e^) are plotted versus the frequency uo and the hyper- 
bolic regions are highlighted in grey. For the chosen fill- 
ing fraction / = 0.4 the first band Ai represents the case 
where Re(e||) < and Re(e^) > whereas in the second 
band A2 the effective permittivities fulfill Re(e||) > and 
Re(e^) < 0. 

To check the pertinence o f the E MT predictions we use 
the exact S-matrix method 



26l-l28j with a finite but large 



number of periods N/2. Beside, for our calculations we 
choose = 100 layers where the last Ge layer is assumed 
to be infinitely large extending to 2; ^ 00. For N ^ 00 
it is well known that propagating waves in a periodic 
structure are Bloch waves and satisfy the Bloch mode 
dispersion relation [26] 
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with Psi = Ps2 = 1, ^pi esic and Pp2 cq, 
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Figure 2: (Color online) Isofreqency curves for the Bloch dis- 
persion relation JS]) (solid lines) and the effective dispersion 
relation Q (dashed lines) for / = 0.4, a — 500 nm and for 
u = 1.55 • lO^^rad/s, 1.6 • lO^^rad/s, 1.7 • lO^^rad/s, and 
1.9 • 10^^ rad/s. 

To analyse the dispersion relation for the exact and 
the effective case we have plotted the isofrequency curves 
in the K-kz plane from Eq. (jl]) and (j5]) for different fre- 
quencies in Fig. [21 Here we used the Bloch Eq. (|5|) 
which gives the exact isofrequency curves for the limit- 
ing case N ^ 00. For the frequencies uj = 1.6-10^^ rad/s 
and UJ = 1.9 • 10^^ rad/s outside the hyperbolic regions 
the effective permittivities are both positive as shown in 
Fig. [TJb) and accordingly we get elliptical isofrequency 
curves. The difference between the principal axes of 
the black and orange curves stems from the fact that 
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Re(e||) < Re(e^) for the frequency range between Ai and 
A2 whereas for frequencies larger (smaher) than A2 (Ai) 
Re(e||) > Re(e^). The Bloch dispersion curves (sohd 
Hne) for these frequency regions fit excellently the effec- 
tive data (dashed line). 

The green and blue curves in Fig. [2] are plotted for 
frequencies which are located in the hyperbolic bands. 
The dashed lines for = 1.55 • lO^^rad/s (blue) and uo = 
1.7 • IQ-'^^rad/s (green) illustrate the hyperbolic effective 
results where the blue curve is an example for a dispersion 
curve from the Ai band and the green one represents 
a dispersion curve from the A2 band. Here the Bloch 
dispersion coincides with the EMT result for small ti. For 
increasing values of n the solid curves show the existence 
of a cj-dependent vertical asymptote as a consequence 
of Bloch Eq. ([5]), while the effective dashed curves do 
not have a maximal ti value. This is a major difference 
between the behavior of the structure predicted by the 
effective theory and by the rigorous theory. 
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Figure 3: (Color online) Log of RC [k < uj/c : 1 — |rp|^; 
K, > (jj/c : Im(rp)] for a SiC/Ge structure and two different k 
regimes with a — 100 nm, / — 0.4 and N — 100. (a), (c) show 
the exact and (b), (d) the effective case. The black solid lines 
indicate the boundaries of the Bloch mode dispersion relation 
from Eq. (|5|). The black dashed lines mark the two hyperbolic 
bands Ai and A2. 

Although the isofrequency curves deduced from the ef- 
fective medium theory and from the Bloch theory coin- 
cide very well for < 1/a as it is well known, the exact 
calculations can deviate very strongly from the effective 
and Bloch results due to the fact that the structure can 
support surface modes which are not taken into account 
neither in the Bloch theory nor in the EMT. To demon- 
strate this fact we have plotted in Fig. [3] both the exact 
and the effective reflection coefficients (RCs) in the uj-k, 
plane in regions of small and large values. In Fig. [3l^a) 
and (c), where the exact RC is plotted, the black solid 
lines represent the boundaries of the Bloch bands given 
by Eq. (j5j) for propagating modes in the medium. Here 



the considered frequency range is still in the first Bloch 
band. Inside these Bloch areas the RC has large val- 
ues and with increasing number N these areas are filled 
due to an increasing number of discrete Bloch modes. 
The maximal value for the contributing Bloch modes 
is therefore given by Eq. ([5]). 

The RC for the EMT is plotted in Fig. [S^b) and (d). 
The black dashed lines mark the hyperbolic bands Ai 
and A2 where Re(e||)Re(e^) < 0. By comparing Fig.j^a) 
with Fig. [3fb) with /^max = 10/a it can be seen that 
the EMT fails to reproduce the exact result for large 
values, as already mentioned before. The boundaries for 
the Bloch bands change drastically whereas the size of 
the hyperbohc bands does not change with increasing k 
values. In Fig. [3jc) and (d) the exact and effective RC 
is plotted up to A^max = Here the RC caclulated 

with the EMT nearly coincides with the exact result. 
The exact calculation shows discrete bloch modes in the 
hyperbolic regime. This comes from the limited number 
of layers {N = 100). For N ^ 00 the discrete Bloch 
modes will fill the whole hyperbolic regime. 

But beyond the Bloch areas in Fig.[3fa) two prominent 
features around uoi = 1.8 • lO^^rad/s (upper branch) and 
UJ2 = 1.6 • lO^^rad/s (lower branch) are visible with very 
high values for the RC. These modes are coupled SPP 
modes stemming from the first SiC layer on the top of 
the structure. These surface modes are damped on a 
scale much smaller than the period a so that they do not 
'feel' the periodic structure. The upper branch around uji 
has a very high RC and moreover high values up to very 
large hi. The SPP contribution becomes therefore rele- 
vant for very small distances z to the periodic structure. 
In the exact calculations, as it can be seen in Fig. [3fa), 
the upper SPP branch is located in the second hyperbolic 
band A2 where the effective results do not show any SPP 
mode in the hyperbolic regimes. The contribution of this 
SPP mode around uji is still visible for < 1/a as it can 
be seen in Fig.[3fc). The contribution of the coupled sur- 
face modes reduces but does not vanish. This qualitative 
difference between both approaches can lead to misinter- 
pretations for very small distances to the structure. 

To estimate the range of validity of the EMT and to 
evaluate the disturbing role played by the SPP mode we 
have finally calculated the local density of states (LDOS) 
of electromagnetic field for all filling fractions at different 
distances z above the SiC / Ge structure both for the ho- 
mogenized structure and for the real one. The results are 
plotted in Fig. 4. For z = 500 nm in Fig.Sl^e) and (f) for 
both plots high LDOS values appear for large / at the 
lower right boundary of A2. This is due to a surface mode 
outside the hyperbolic bands where Re(e||)Re(e^) = 1 
and Re(e||) < 0. These conditions together with Eq. (|2]) 
lead to /> 0.5 as a necessary condition for having sur- 
face modes outside the hyperbolic band within the EMT. 
The contribution of the coupled SPP mode of the first 
layer within the hyperbolic band A2 is very low and the 
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Figure 4: (Color online) Local Density of States (LDOS) in 
the Lu-f plane for different distances z above the SiC/Ge bi- 
layer structure calculated with the exact (a), (c), (e) and the 
effective (b), (d), (f) theory with a = 100 nm and N = 100. 
The dot-dashed white lines mark the hyperbolic regions Ai 
and A2. 



effective calculations fit very well the exact results. How- 
ever, when decreasing the distance z this contribution of 
the coupled SPP mode becomes larger and eventually for 
z = 50 nm the effective LDOS is not able to mimic the 
exact LDOS as illustrated in Fig. HJ^a) and (b). Here the 
contribution of the coupled SPP mode of the first layer 
is much larger than the contribution of the Bloch waves. 
Hence, from the EMT calculations one would expect a 
large LDOS stemming from the hyperbolic modes, which 
would be a wrong conclusion, since as we see from the 
exact results the large LDOS is due to SPP modes from 
the first layer. 

In this work we have shown that even for distances 
larger than the lattice period the EMT does not provide 
a correct physical picture of the structure when localized 
modes such as surface polaritons are present. We have 
shown that the usual conditions for the validity of the 
EMT a <C A and /^a <C 1 have to be augmented by the 
condition that the distance z at which the electromagntic 
field is investigated must be larger than the penetration 



length of the SPP modes into the vacuum region. 
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